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ABSTRACT 

For a wide class of self-gravitating systems, we show that if the density is cusped like r near the center, then 
the limiting value of the anisotropy parameter /3 - I - (v^)/(2(v^)) at the center may not be greater than y/2. 
Here (v^) and (v^) are the radial and tangential velocity second moments. This follows from the nonnegativity 
of the phase-space density. We compare this theorem to other proposed relations between the cusp slope and 
the central anisotropy to clarify their applicabilities and underlying assumptions. The extension of this theorem 
to tracer populations in an externally imposed potential is also derived. In particular, for stars moving in the 
vicinity of a central black hole, this reduces to j > /3 + j, indicating that an isotropic system in Keplerian 
potential should be cusped at least as steep as r Similar limits have been noticed before for specific forms 
of the distribution function, but here we establish this as a general result. 

Subject headings: stellar dynamics — methods: analytical — galaxies: kinematics and dynamics 



1. INTRODUCTION 

Numerical simulations of halo formation provide strong ev- 
idence that the inner parts of dark matter halos are strongly 
cusped. Typically, the density profile p behaves like_r_^ 
where y lies between 1 and 1.5 (Navarro. Frenk. & White! 
119951 iMoore et al."1998V Although this numerical result 
seems well established, observational evidence that dark ha- 
los are cusped has been surprisingly elusive. A disparate 
body of data - including the rotation curves of dwarf spi- 
ral galaxies (Palunas & WilHams 2000; de Blok et al. 2001), 
the kinematics of Local Group dwarf spheroidal galaxies 
(jfClevna et al. 2003), and mass models of gra vitational lens 
systems (Tvson. Ko chanski. & deirAntoniolll998ll) - seem to 
favor constant-density cores. 

Recently, Hansen (2004) claimed that the only density 
slopes permitted by the spherically symmetric and isotropic 
Jeans equations are 1 < y < 3 if the phase-s pace density-hke 
quantity, p/{v^)^^^, follows a power law ( Tavlor & Navarrd 
12001). This result was inferred from the condition that 
the power-law solution of the Jeans and Poisson equations 
is physical, subject to the "equation of state" (EOS), p oc 
r"''(v^)^, where p and c are constants. He further argued 
that, for the system with an anisotropic velocity dispersion 
tensor, 1 -H yS < y < 3. Here the anisotropy parameter is 
P - I - {Vj)/{2{v^.)), where (vj) and (V j) are the radial and 
the tangential velocity second moments llBinnev & Tremain^ 
[1287). Given the nearly isotropic conditions found in the cen- 
tral parts of simulated dark halos, this already seems to indi- 
cate that the density profile cannot be shallower than p ~ r ' . 

The idea of looking for constraints between the central den- 
sity slope and the anisotropy is an excellent one. In § |2l 
we show that the inequality y > 2/? is a necessary condi- 
tion for the nonnegativity of the distribution function (DF). 
This generalizes two well-established results: (1) a spherical 
system with a hole in the center cannot be supported by an 
isotropic velocity dispersion tensor (i.e., y > 2y6 = 0; see 
[iremaine et al. 1994), and (2) a spherical system with parti- 
cles in purely radial orbits cannot be supported by a density 
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cusp shallower than the isothermal cusp (i.e., y > 2/3 - 2; see 
Richstone & Tremaine 1984). In §|3l w e stu dv the scale-free 
power-law cusps that inspired Hansen (2004) and show that 
the inequality y > 1 +/? is related to the boundary condition at 
infinity rather than at the center. Finally, in § |3 we consider 
the generalization of the theorem to tracer populations in an 
externally imposed cusped potential - for example, to stars 
moving around a central black hole. In particular, with the 
potential cusped as r (0 < 6 < 1), we find that the inequal- 
ity becomes y > 2/3 + - (3)5. While a similar limit has been 
noted in the literature for the case in which the DF is scale-free 
(White 1981; de Bruiine, van der Mai-el, & de Zeeuw 1996), 
we derive the limit as a general condition for a DF with a cen- 
trally diverging (and not necessarily self-consistent) potential. 

2. PROOF OF A CUSP SLOPE-ANISOTROPY THEOREM 

Here we shall prove that if the density is cusped like r 
near the center, then the limiting value of /3 at the center may 
not be greater than y/2. This relation is indeed suggested by 
solving the Jeans equation for constant /3. That is, the one- 
dimensional radial velocity dispersion obtained as the solu- 
tion, in general, diverges at the center if y < 2/3, which is un- 
physical unless the central potential well depth is infinite. In 
fact, we will show that the inequality y >2/3\s the necessary 
condition for the nonnegativity of the DF. 



2.1. Constant Anisotropy Distribution Functions 
First, let us suppose that the DF is given by the Ansatz^ 

f{E,L)^L-^PfE{E), 



(1) 



where L is the specific angular momentum, and fsiE) is a 
function of the binding energy E alone. This seems as though 
it is a restrictive assumption, but this is not really the case. 
Rather, equation ([3 arises naturally as the most divergent 
term in a Laurent series expansion with respect to L at L = 
for a very wide class of DFs. By integrating equation (Q over 

^ Strictly speaking, the Ansatz is valid only for /} < 1. However, it can be 
extended to /J = 1 using the relation lim/j^, [L^^^/ni -/3)] = S(0, where 
(5(.y) is the Dirac delta "function." Subsequently, equations l2l isl are still 
valid for f(E, L) = S(l?)fE(E) as thcyS = 1 limit simply without the r(l 
factor 
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velocity space, we find that tlie density is given by 

p = r-^^Dp fV - E)'l^-^fE{E) dE, 
Jo 



(2) 



(3) 



where if/ is the relative potential and r(x) is the gamma func- 
tion. The unknown function fsiE) can be fou nd from the for- 
mula (Cuddeford 1991; Kochanek 1996; Wi lkinson & EvansI 
fr999.;.An & Evans.2006: Eivans & Aa2006.) 



fE(E) = Cp 



Jo 



dxli d"+^h 1 d"h 
E" dtf/" 



(E - tf/y #"+1 
^ 

(27r)3/2r(l -a)r(l -jS)" 



(4) 



(5) 



Here h = r^p is expressed as a function of i/r, and n - [§ -y6J 

and a - |-j0-nare the integer floor and the fractional part 

of I - /?• It is a simple exercise to show that the anisotropy 
parameter for this model is the same as the parameter /3 in the 
expression of DF (eq.[0. 

By considering the limit of equation Q as r ^ one can 
infer that p should diverge at least as fast as if /5 > 
and cannot approach zero faster than r^'^'^' if y6 < unless 
the integral vanishes in the same limit. This in fact is the 
basic argument that leads to the theorem. In the following, we 
provide a more stringent proof of the theorem. 

2.1.1. The Case/3< \ 

For 1 3 < 5, a direct generalization of the proof of 
[Iremaine et al. ( 1994) for isotropic models suflices. That is, 
from equation (0, we find that 



d - rjEiEUE 



2/*r( 1/2-/3) 



(6) 



(7) 



for any physical DF, as the integrand is always positive. Sim- 
ilarly, if /? = 5, then 

rfE(E)dE^2jT^fEW>0. (8) 
Jo 



— (rp) = 2n^ — 
dijj dij/ 



However, 



dh 



dr hd\nh 



Pl2f,^dlnp 



(9) 



dip r din r difi/dr\ dlnr 
Since dtfr/dr < for any physical potential, we thus find that 

dlnp 



2p<- 



d\\\r 



(10) 



This holds everywhere. Specializing to the limit at the center, 
we obtain the desired result. 

2.1.2. The Case/3> \ 

When P > \, equation (|6j is invalid, and therefore a differ- 
ent proof is required. For this purpose, we first note that 



>0 



(11) 



for < E < i// and (3 > \- Then from equation (|2|i, we find 
that 



P> 



Dp r 



fE{E)dE. 



(12) 



Therefore, 



lim (pr^V"''^) j ME) dE>0, (13) 

where iffQ = iff{r - 0). Next, if we consider the case in which 
tf/Q is finite (see §|4]for the centrally diverging potential), we 
have 

limfpr^V"'^') = i/t"^lim(pr^'') > 0. (14) 

r— >0 ^ ' r— >0 ^ ^ 

Since 1^0 > 0, we find that lim,.^() (pr^''') > 0, that is, lim,.^() h 
either is nonzero positive and finite or diverges to positive in- 
finity. If the former is the case, it is straightforward to show 
that (using FHopital's rule) 



hr d(hr) d ln(hr) 

lim h — lim — = lim — ; — = lim h - 



r^O r^O r r^O dr 

and consequently that 

d\n{hr) 

lim = 1 

r^o dlnr 



^0 dlnr 



,. dlnh 
lim — — = 0. 



^0 dlnr 



(15) 



(16) 



On the other hand, if limr^oh is divergent, I'Hopital's rule 
indicates that 

dlnh Inh 
hm— — =lim— . (17) 
r^o d In r r^o In r 

However, we have (Inh/ Inr) < for sufficiently small r, so 



In/i 

hm < 

r^o In r 

Hence, for both cases. 



dlnh 
lim — — < 0. 
r^o li In r 



dlnh dlnp 
lim — — = lim — +2B<0 
i-^Q dlnr r^o dlnr 

,. dlnp 
7 = - hm > 2/3, 



^0 dlnr 



(18) 

(19) 
(20) 



which is the desired result. 



2.2. General (Analytic) Distribution Functions 

This result is in fact far more general than the assumed 
Ansatz for the DF (eq. 0. For example, Ciotti & Pelleariml 
(1992) fou nd the same l imit (y > 2/3 - 0) for the Osicko^ 
(ll979a,b.) - Merritll (I1985h type DF, which is manifestly not in 
the form of equation ([T). We extend the limit derived in the 
preceding section to a wide class of DFs by the following sim- 
ple argument. In general, any analytic DF can be expressed 
either as sums of equation Q or in terms of a Laurent series 
expansion with respect to L at L = (really a special class of 
the former). Then since L = rvj, as r —> Q, the DF is domi- 
nated by the term associated with the leading order of L, and 
consequently so is the behavior of the density near the center. 
It is also straightforward to show that the anisotropy parame- 
ter at the center is indeed determined by /3o, where '-2/3{) is 
the power to the leading term of L. The desired result there- 
fore follows from the preceding proof for the special case of 
the DF with a single term. In Appendix, we discuss the con- 
ditions of applicability of the proof more mathematically and 
argue that the theorem holds for all physically reasonable DFs 
of spherical systems. 

We find that the cusp slope-anisotropy theorem (y > 2/3) is 
actually quite reasonable. It implies that, if the anisotropy is 
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radially biased (J3 > 0) near the center, the density is cusped, 
and that, unless the cusp slope is steeper than that of isother- 
mal cusp (y = 2), there is a finite upper limit to /3 that is 
strictly smaller than unity. Similarly, if the density is cored 
(y - 0), the central anisotropy is either tangentially biased or 
at most isotropic. 

3. SCALE-FREE DENSITY PROFILE 

Recently, iHansenl i2004l) derived a similar but stricter in- 
equality y > 1 + /3, based on the condition that there exists 
a physical power-law solution to the spherically symmetric 
Jeans and Poisson equations with the EOS, p oc r^''{v^^y. 
Since /3 < 1, his result is stronger than our result. However, 
we note that his result is strictly valid only if both p and (v^) 
behave as the pure power law extending globally to the infin- 
ity. In fact, we find that the condition y > /3+1 originates from 
the boundary condition at infinity rather than at the center and 
thus argue that the result should be understood to refer to the 
asymptotic density power index at infinity, not the central den- 
sity slope. In particular, the supposed piece of evidence cited 
by Hansen (2004) for his inequ ality (Lokas & Hoffman 2003 
see also B ouvier & .T anin 1968) is in fact due to the constraint 
imposed by y > 2yS on the central slope through the positivity 
of the DF (§|2j see also Richstone & Tremaine 1984), since 
it involves the case of purely radial motion (/?=!) for which 
the two limits coincide {y > 2). 

The general integral solution of the Jeans equation with 
constant /? can be written by admitting an integration constant 
B, 

r2Vv?) = B+ [' dr' r'^^pir') ^ . 

Jrr, ar r=r' 



(21) 



The potential gradient (i.e., the local gravitational accelera- 
tion) can be found from the enclosed mass 

GMr 



dr 



•Jrn 



Mr^Mr,+An\ p{r')r"-dr', (22) 



where the negative sign is due to our choice of sign for i/r. 
If we assume a strict power-law behavior for the density as 
p - Ar^^, where A > is constant, then is given by (y + 3) 

M, = M,„ + 47rA(3 - yY' [r"-^ - r]-'') . (23) 

If the power law extends to the center, the result must be valid 
for the choice of ro = 0. However, then the mass within any 
finite radius always diverges for y > 3 even if Mq is finite, 
and therefore the model is unphysical. The resulting upper 
limit y < 3 is well established. By substituting equation \22\ 
and p - Ar^'*' in equation (12 1> . (v^) as a function of r (y 9^ 3, 
y + 13+ I, and y 2/3 - 1) is found to be 



(3 - y)(fi ■ 



1 

+ - 
r 



47TGAr: 



y+1) 
GM,, 



(3-y)(2yS-y-l) 2^ - y - 1 



(24) 



where B is an integration constant to be determined from the 
boundary condition. Here if we assume strictly scale-free be- 
havior for p, equation (I24> is valid everywhere extending from 
r = to r = oo. With ro = and Mq = 0, we can show that the 
condition for (v^) to be nonne gative everywh ere leads to the 
inequality y > as found by Hansen (7004), and Z? > 0. In 
addition, the self-similarity implies strict power-law behavior 
for (v^) as well. Since A > 0, this can only be obtained with 
the choice of B = 0. In fact, the choice of B = can inde- 
pendently be deduced from the boundary condition at infinity. 



That is, (v^) is nondivergent for a finite potential (y < 2), or if 
the potential diverges (y > 2), the velocity dispersion cannot 
diverge faster than the potential does. 

However, if we relax the assumption that p is strictly scale- 
free everywhere and replace it with p being locally well ap- 
proximated by the power law near the center, equation (I24t 
is valid only for the region where p ^ Ar^~^ , so the condi- 
tion that (v^) is nonnegative only needs to be checked for this 
regime. Provided that the power law provides a good approx- 
imation to the behavior of p near the center and that we limit 
our attention to a self-consistent density-potential system, it 
is reasonable to choose ro = and Mo - 0. Then for y > 2yS, 
with any positive constant B, equation (I24> returns the correct 
behavior of the velocity dispersion near the center, although 
its validity does not extend to infinity. For y > y6 H- 1, the 
velocity dispersion (ro = Mo = 0) is given by 



<v?> 



27rGA 



(3-y)(y-/?-l) 



(25) 



Near the center, we find that (v^) ~ A'r^"'', where A' is a 
positive constant. This is just the approximate local power- 
law behavior of the velocity dispersion near the center (valid 
locally for any B). In addition, if y6 -i- 1 < y < 2, the cen- 
tral velocity dispersion is finite. Although (v^) diverges at the 
center ify>2>y6+l, this behavior can be physical be- 
cause the self-consistent potential well depth for this case is 
also infinite. On the other hand, if y < y6 + 1, the same velocity 
dispersion can be written as 



B _ ^203+ i-r). 



2;rGA 



(3 - y)(j8 + 1 - y) 



(26) 



In other words, since y-2yS<2-yify<j6-Hl, the lead- 
ing term for the velocity dispersion near the center is given 
by <v^> ~ Bf'^'^f^ provided that B 0. While not strictly 
scale-free, the local behavior of the velocity dispersion near 
the center can still be approximated as a power law, and fur- 
thermore, as long as B > and < y - 2/3 < 2 - y, it is 
positive and finite. While we note that for a sufficiently large 
r, equation ( I26t eventually becomes negative, since the r^"'' 
term becomes dominant as r ^ oo, this does not restrict the 
central power index for density, provided that the behavior of 
p starts to deviate from ~ Ar^^ (toward the steeper falloff) at 
smaller r than the value at which (v^) = in equation ( [26t . 

Let us consider the explicit example of the Hernquist ( 199Cj) 
model, which has a r"^ density cusp (i.e., y = 1). The ra- 
dial velocity dispersion of the constant-j6 m odel is (see, e.g., 
iBaes & Deionghel2002tlEvans & AiJ200l 



GM 



(5 - 2/3) (r -H fl)2 



2^1 (l, 5; 6 -2/?;^], (27) 
\ r + al 



where the potential is given by = GM/(r + a). It is straight- 
forward to show that, for Q < p < i, (v^) is everywhere posi- 
tive finite with its leading order behavior near the center given 



by 



.1-2/3 



(ifyS 



1 



2, then (v^) = (A/4; see also IE vans 
I2005h . On the other hand, if y8 < 0, the leading order for <y^) 
near the center is found to be ~ r with a positive coefficient. 
For the isotropic case {fi - 0), the elementary functional ex- 
pression for t he vel ocity dispersion is given in equation (10) 
of Hernauisll ( 119901) . whose lead ing order beha vior is found to 
be ~ rlnr"' (see, e.g., eq. 11 of Hernquist 1990). 

As another example, we consider the dark matter profile 
proposed by De hnen & McLau ghlin (2005). They solved the 
spherically sy mmetric Jeans a nd Poisson equations with the 
same EOS as in lHanseJ ( 12004ft . However, they found a family 
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of "realistic" models with a finite mass and an infinite extent. 
With p oc r^P{v^y^^, the inner density cusp of their models 
is given by y = (7 + 10y6())/9, where 0q is the anisotropy pa- 
rameter at the center. It is clear that for all members of these 
models y < y6o + 1, since ySo < 1 < 2, thus violating the in- 
equality y > /3 + I. On the other hand, our result y > ip 
indicates that they are physical if and only if Pq < ^. \n fact, 

(v^) near the center for this family behaves as ~ r*^^"^^"*^', and 
thus the limit y > 2/? is equivalent to the condition that the 
central velocity dispersion is finite. 

A similar analysis of equation (I24> can be applied to dis- 
cover the asymptotic behavior of the velocity dispersion at 
infinity. Suppose that p asymptotically approaches a power 
law for a sufficiently large r. Then for the same range, equa- 
tion ( I24> is a valid expression for (v^), provided that the power 
law ~ Ap'~>' describes the asymptotic behavior of p and M^j, is 
the mass within tq. Here B - Q from the boundary condition 
at infinity. If the total mass is finite {y > 3> > fi+ \ > 2j3 - 1), 
we can simply set ro = oo and Moo - M, where M is the total 
mass. Then since -1 > 2 - y, the asymptotic behavior of the 
velocity dispersion is given by (v^) ~ (y -i- 1 - ipY^GM j r, 
which is consistent with Keplerian falloff. For an infinite- 
mass system, we find that 2 - y > -1, and thus the leading 
term of (v^) for r ^ oo is ~ A'p-^^ . Here A' is a positive con- 
stant if y > yS H- 1 , whereas it is a negative constant if y < j6 -i- 1 . 
In other words, from the condition of nonnegativity of the ve- 
locity dispersion toward infinity, we actually recover the in- 
equality of Hansen (2004), y > /? + 1, although here y and fi 
are the asymptotic density power index and the anisotropy pa- 
rameter at infinity. Although the velocity dispersion diverges 
as r ^ oo if y < 2, the behavior may be acceptable because 
the potential also diverges with the same power index, so the 
system is still bounded. 

4. CENTRAL BLACK HOLE 

We have found one further exampl e of a cusp slope- 
anisotropy relationship in the literature. White (1981) found 
a relation** by studying the special case of scale-free densities 
in scale-free potentials, namely 



y>-+P(2-6), 



(28) 



where 6 is the central power-law index for the potential (i.e., 
ill ~ which ma y be externall y imposed. Note the change 
of the notation from Whi't^ (1198 Ih . The form of the limit given 
in equation ( I28t is actually that of de Bruiine et al. ( 199^, 
who performed a similar study to White (1981) but allowed 
for flattening. This result was derived from a specific func- 
tional form for the DF, and in particular assumed that fE{E) is 
scale-free. We note that, since yS < 1, the inequality y > 2yS is 
automatically satisfied if y > 2. For y < 2, the self-consistent 
potential-density has a finite central potential, so (5 = 0. In 
this case, equation (I28> reduces to y > 2yS. However, in the 
presence of a central black hole {6 - 1), equation ( I28> pro- 
vides us with a different limit y > /? + 5, which is stricter than 

y>2ySify3< i. 

Here we note that this limit can, without assuming that the 
density or DF is scale-free, be derived from the nonnegativity 
of the DF for a massless tracer population in the Keplerian po- 
tential of a central point mass. Let us suppose that the DF for 



The inequality given in'Wliits? (1981") or'de Bruiine et alj 11993) did not 
include the case of equality because of the specific form of the scale-free DF. 
The general result actually extends to include the case of equality through the 
transition of the power-law distribution to the Dirac delta distribution. 



a massless tracer population is given by equation (^. (The 
results can then be extended to more general DFs using the 
identical argument of § 12.21 ) Next, we assume that these trac- 
ers are subject to the potential = GMj r of a point mass at 
the center. Then the number density of the tracer population 
n can be found from integration of the DF over velocity space 
as in equation ^\ 



"";^Jo \—~^} f^^^^'^^ 
^.1/2 j„ 



-GMIr 



(GM -rEy^^-PfE{E)dE 



For y6 < 5, we find that 



, pCM/r 

dr^ > "Jo ( 



EfE(E)dE 



(GM - rEfl^+P 
for any nonnegative DF. Consequently, 

or, equivalently. 



d\nn 
liln r 



<0 



< 



d\nn 1 

>B+- 

d\nr ^ 2 



(29) 



(30) 



(31) 



(32) 



which is the desired result. In fact, if the external potential 
is replaced by i/f - Cjr^ in equation M9\ . it is easy to see 
that equation ( 12 8> simply follows from an essentially identical 
argument with y - -(dlnn/dlnr) being the central power 
index for the number density of the tracer population. 

In fact, the result can still be obtained without assuming the 
specific form of (scale-free) potential. That is, provided only 
that 6 - - lim,-^Q{d Inilf / d In r) and < 6 < 1, we find from 
equation (|5J that 



D, 



-A 



1/2-/} nip 



„2/3 



Jo 



1-^ 



fE(E)dE, 



(33) 



and thus for y6 < 5, we have 
(nr 



Jo ' 



d 

difr ' 

However, of course. 



d I 
— [nr 



-1/2 



EfE{E)dE 

„^2/J-l^/J-I/2 ^ 



> 0. (34) 



3/2 



d\n t///dlnr 
and by taking the limit r 



di///dr dlnr 
dlnn I 1 

d^^r^'^^^-2 



In (n 

dlnij/ 
dlnr 



-l/2\ 



> 0, (35) 



0, 



r>2y8 + (- 



/3 S. 



(36) 



Here note that linir^oid In iff /dlnr) = -6 < Q. This also indi- 
cates that the result is still valid even if the self-gravity of the 
tracer population is appreciable, as long as the potential (that 
may be decoupled from the density) is divergent at the center 
If j6 = J, the limit given in equation ( I28l l is, in fact, identical 
to y > 2yS = 1. Since the derivation of the limit in §|2]for /? - 
5 does not use the assumption of self-consistency, it is still 
applicable here. Therefore, the limits in equations ( I32> and 



36) can be extended toyS < 4. For £< ^ < 1, the limit y > 2/? 



is actually stronger than equation (I28> (note that < 6 < 1) 
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However, with a ce ntrally diverging potential (i.e., i/'o = oo), 
the proof given in § l2.L2l is not directly applicable. Instead, 
from equation il3\ . we now find that linir^o^^ where k - 
nr^i//^^/^ either is finite and nonzero or diverges to p ositive 
infinity. Following exactly the same argument as in § I2.L2I 
applied to k = nr^i//^^^^ instead of h - pr^^, we have 



lim 



dink 



^0 dlnr 



d\nn I 1 

lim— — +2yS + r6- - 



lilni^ 
lim — — < 0, 
r^o d\nr 



(37) 



which translates to equation ( I36t . and thus we can extend the 
limit of equation \36\ to /3 > ^ as well. 

The limit y > /3 + j for 5=1 indicates that a spherical 
isotropic system subject to a Keplerian potential should pos- 
sess a central density cusp at least as steep as r"'^^. Similarly, 
if an isotropic stellar system is subject to a divergent dark mat- 
ter potential due to a cusped profile with a slope steeper than 
that of the isothermal cusp, the stellar system should also be 
cusped with its cusp slope constrained to be y*. > (701^/2)- 1, 
where is the cusp index for the stellar system and yDM > 2 
is that of the dark matter profile. Of course, if y^M < 2, the 
central potential is finite provided that there is no other source 
of divergent potential, and thus the limit simply reduces to 
y* > 0. On the other hand, if the system were mildly radially 
anisotropic (near the center), that is to say /3 ~ j, the limit for 
the supportable cusp slope would be steeper, much like y > I. 

5. CONCLUSIONS 

We have proved, for a very wide class of steady-state gravi- 
tating system, a theorem constraining the central cusp slope 



of the density profile y (eq. I20> and the central velocity 
anisotropy Specifically, the inequality y > 2y6 is a neces- 
sary condition for the nonnegativity of the distribution func- 
tion (DP). If there is a divergent external potential, decou- 
pled from the density profile, then the inequality generalizes 
to y > 2/3 - (J3 - ^)6. Here the external potential diverges as 
~ r at the center. Finally, if the external potential is due to 
a central black hole, it reduces to y > yS + i. We expect our 
result to be useful in the study of dense stellar systems, or in 
the building of extreme stellar dynamical models. 

As most N-body simulations predict only very modest 
anisotropics (J3 ^ 0) in the very center, the application of 
our result does not directly constrain the central density pro- 
file (y > 0). While the inequality derived by Hansen (2004), 
namely, I + /3 < y < 3, appears to be stronger than our re- 
sult, his lower limit is only strictly applicable to the scale-free 
power-law density profile of infinite extent. It appears that 
y > /3 + 1 actually constrains the asymptotic behavior of the 
density power index and the anisotropy parameter at infinity 
rather than at the center 
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APPENDIX 

Here we provide a more detailed argument for the generality of the theorem than that given in § 12.21 While we do not claim that 
the following argument strictly adheres to the high standard of the pure mathematician, we hope that it indicates the generality of 
the result. 

First, let us suppose that a DF f(E, L) can be written as 



f{E,L)^L-^l'-{UE) + ,m,L)], 
where /o(£') is a function of E alone, whereas f\(E, L) is a continuous function that satisfies 

/i(£,L = 0) = 0, 

which further implies that /o(£) > for all accessible values of E from the nonnegativity of the DF. Then 

P^^^jj l^'^"' ihiE) + /i (E, L)] VjdvjdVr 

p{vj}^27T JJ vjL-^f" [fo(E)+ME,L)]vTdvTdv, 
p<v4 > = 27r ^ v4l-2A [/o(£) + /i {E, L)] vjdvjdv. 



(1) 
(2) 



4D^„(1-A)) 
■r2A(3-2^o) 



jyEi^-Ef'^-^^ME) + ^ ff dvjdv,v\ 



■0 f 
Ji 



v\ + vl 



(3) 



(4) 



(5) 



Taking the limit r — > 0, the velocity moment integrals of the /i term vanishes, provided that the domain over which it is integrated 
is bounded and that /i is sufficiently well behaved. Note then that the anisotropy parameter at the center for this model is found 
to be ^ 

lim^=l-lim^^ = l-(l-,6o)=A)- (6) 



6 AN & EVANS 

Next, suppose that ho = r^^''po and po is the density profile built by the DF, f(E, L) - L^^''fo(E). Then the proof given in § 12.11 
indicates that yo > 2y6o, where 

ro = -lim— . (7) 

r^o dlnr 

Now, since lim,.^()(t/ln/io/iilnr) = 2/3o - yo < 0, we have linir^o/io 9^ 0, that is, lim^^o/io is nonzero finite or diverges. 
Subsequently, from equation 



^-1/10 = 0, (8) 

ho J 



but lim,.^o ho + 0, and therefore we find that 



lim— ^ = 1. (9) 



If limr^o ho is finite (i.e. yo - 2y3o), then equation (|9} implies that 

< limr^% = lim/io < oo, (10) 

r^O r^O 

SO lim,.^o '"^'^"p is also finite, that is, p ~ r"^'^". On the other hand, if limr^o ho - oo, I'Hopital's rule indicates that 
. r^^V d(r^<'p)ldr r'^'^" p d\n(r^° p) I d\n r 1 / d\np\ 

lim — — - hm — - — — — = lim — — ; — — = 2Bo - hm — - — = 1 ■ (11) 

r->o «o '■^0 dho/dr r^o ho dinho/dinr 2po-yo\ r^o d In rj 

In other words, combining the results in equations (I10> and (O, we have established that 

,. d\np d\npo 

hm— — =yo = lim— I , (12) 

r^o dinr r^o dinr 

where p is the density profile built by the DF of equation ([0 and po by f(E, L) - L^^^^ foiE). Finally, taking equations and 
J12> together, we have established that the theorem yo > 2Po extends to a more general class of DFs of the form of equation ([Q. 

How general is the form of the DF in equation O? We argue that it is almost always possible to reduce most well-behaving 
DFs to the form of equation ([Q. That is, for a general DF of a spherically symmetric system in equilibrium, the reduction is 
possible if there exists a <2 such that 

fE(E)^\imL"f(E,L), (13) 

where fE{E) should be finite and nonzero for values of E at least in some nonempty subset of all the accessible values of 
E e [0, i/'o]. Then the original DF can be written in a form of equation Q as 

/(£, L) = L " IfE(E) + ME, D] ■ ME, L) = Vf{E, L) - fsiE). (14) 

and it is obvious to show that 

/i {E, L = 0) = lim L"f{E, L) - fsiE) = 0. (15) 
For example, for the DF of the form given bv lCuddefordldTOQlh . 

/(£, L) = L-2/5»/o(0 ; Q = E-^, (16) 

2ri 

the reduction is given by 

f(E, L) = [ME) + /i (E, D] ; /i (E, L) = MQ) ' ME). (17) 
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